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Abstract 



A massless Dirac particle is considered, moving along the x-axis 
while Pauli-coupled by its anomalous magnetic moment to a piece- 
wise constant magnetic field along the same axis, with stochastically 
varying sign. The motion is approximated as a quantum walk with 
unitary noise, for which the evolution can be found exactly. Initially 
ballistic, the motion approaches a classical diffusion on a time-scale 
determined by the speed of light, the size of the magnetic moment, 
the strength of the field and the time interval between changes in its 
direction. It is suggested that a process of this type could occur in 
the Sun's corona, significantly affecting the solar fluxes of one or more 
neutrino types. 
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1 Introduction 



It is known that the time-evolution of some free-particle relativistic wave 
equations can be approximated by a quantum walk (QW). This connection 
has its origins in Feynman's path-integral approach to the propagator for 
Dirac's equation |U |2J, and so predates the extensive researches into QWs 
and their potential for applications in quantum information theory that have 
followed from the seminal paper of Aharonov et al. [3]. 

The connection has been explored in work by Childs and Goldstone on the 
massless Dirac equation |4J, by Katori et al. [5] on Weyl's neutrino equation, 
and by Strauch [6] and us [7] on Dirac's equation with a non-zero mass. As 
a result it is now well-understood how the dynamics of a free relativistic 
particle with spin 1/2, moving on a line, can be approximated arbitrarily 
closely by a simple, one-dimensional QW. 

Because the free Weyl and Dirac equations can be solved easily and ex- 
actly, the connections with QWs are mathematically interesting but of lim- 
ited importance to relativistic physics. In contrast, we consider here a rela- 
tivistic quantum system that does not seem to be amenable to analysis unless 
approximated as a noisy QW. 

Hackett [8] considered the effect of adding an arbitrarily small amount 
of classical randomness in a unitary way to a simple QW on the line. The 
idea of a QW contaminated by unitary noise was then explored extensively 
by Shapira et al. |9J. An initially ballistic motion with a spreading rate 
proportional to the elapsed time, as typical of a QW, is eventually replaced 
by a diffusive motion, with a spreading rate proportional to the square-root 
of the time, as typical of a classical random walk (CRW). The transition 
occurs after a time (number of iterations) determined by the strength of the 
noise. Addition of noise in other forms to QWs has been shown to produce 
similar effects [TQ1 [EE] . 

Recently we have studied a particularly simple one-dimensional QW with 
unitary noise, with the property that an initially diagonal density matrix 
remains diagonal during the evolution of the system [121 03] • The form of 
this evolving density matrix has been found exactly. This simple system may 
be regarded as a toy model for QWs with unitary noise, and it does exhibit 
the transition from ballistic to diffusional behavior explicitly. 

What we now show is that, despite its simplicity, this model has an appli- 
cation to the description of a neutral, massless Dirac particle with anomalous 
magnetic moment, Pauli-coupled to a piece- wise constant magnetic field with 
stochastically changing direction, and confined to move on a straight line. 

The interpretation of the transitional behavior of the QW in the context 
of this relativistic system is rather remarkable: at short times t there is a 
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high probability of finding the particle at a distance ct from its starting point, 
where c is the speed of light, as expected for a free massless particle. But 
as time passes, the motion tends more and more towards a diffusion on the 
line, with a high probability of finding the particle within a distance a\/t 
of a certain point xq as a normal distribution is approached, centered on 
that point. Here a is a constant whose value, together with that of xo, is 
determined by the speed of light, the magnetic moment, the strength of the 
magnetic field, and the time interval between changes in field direction. 

2 The relativistic system 

The dynamics of the particle in this case is governed by Dirac's equation in 
the form 

ihdip(x,t)/dt = H(t)ip(x,t) , 

where H(t) = cd ■ p- fi S ■ B(t) . (1) 

Here p = —ihd/dx is the momentum 3- vector, \x is the magnetic moment, 
B(t) is the external magnetic field 3- vector, and S is the spin 3- vector. The 
4x4 matrices in ([!]) are conveniently defined in terms of the 2x2 Pauli 
matrices a and 2x2 unit matrix I 2 by 

a = 0-3 ®a, (3 = a 2 ®I 2 , S=\hl 2 ®a. (2) 

Can such a system be realized in Nature or the laboratory? We postpone 
a discussion of this until the end of the paper; at this stage our objective 
is to add to the collection of relativistic systems amenable to mathematical 
analysis. 

To proceed, suppose that p 2 \tp) = Pz \i>) = on states \ij)(t)) = ip(x,t) 
of interest, writing x\ = x, pi = p = —ihd/dx. We consider magnetic fields 
of the form B(t) = (B(t), 0, 0) with 

B(t) = sB , NA < t < (N + 1)A, iV = 0, 1, (3) 

where A > and B Q are constants, and s — ±1, each sign having probability 
1/2, independently at each value of N. Under these assumptions, the helicity 
operator |/2®o"i is a constant of motion, and we may suppose for definiteness 
that 

(iaOo-OlV) = |^>, (4) 
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and work henceforth in this eigenspace, neglecting the action of the Pauli 
matrices in the second factor of the tensor product. Now the Hamiltonian 
reduces effectively to 

H = cpa 3 -(hfMB(t)/2)<T 2 , (5) 

and the evolution operator carrying \ip(t = nA)) into \ip(t — (n + 1)A)} 
reduces to 

V{s) = exp{-i{cpa 3 -s(hfiB /2)a 2 }A/h}. (6) 



3 Approximation as a noisy quantum walk 

The key step in our approximate treatment is to write from ^ 

V(s) « exp{-i[cpA/h]a 3 } exp{i[s^B A/2]a 2 } • (7) 

Bearing in mind the Campbell-Baker-Hausdorff formula [2] , we assume that 
([7]) is a good approximation provided 

\c(p) n A/h\ < 1 and \tiB A/2\ < 1 , (8) 

where {p) n is the expectation value of p in the state \ip(t = nA)). We shall not 
attempt here a more rigorous analysis of the approximation §7§ , analogous to 
that given for the free-particle [7] , but content ourselves with the assumption 
that inequalities ^ hold so strongly for each n = , 1 , 2 , . . . ,N that there is 
a negligible accumulation of errors when the approximate evolution operators 
in the form ^ are applied N times to \ip(0)), while the system evolves over 
a time t = NA of interest (see below). 

The operators (j5J), ^ and ^ act on the Hilbert space H = Hw ® He, 
where Tiw is the space of square- integrable functions of x on which p acts, and 
He is the 2-dimensional complex vector space on which the Pauli matrices 
in ((5]) and ^ act, each space having the usual scalar product. In Hw we 
introduce the infinite sequence of orthonormal states \k) defined by 

|Jfe) = F(a;-JfecA), k = , ±1 , ±2 , . . . , (9) 

where F(x) is some chosen smooth function with compact support S < 
(-cA/2,cA/2), satisfying 

J \F(x)\ 2 dx = l. (10) 
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Note that the states \k) do not in general form a basis in Hw We also 
introduce on Tiw, the translation operators 

E± = exp (TicpA/h) , E±\k) = \k±l). (11) 

In He we introduce the orthonormal basis of states |r), t = ±1, with 

<7 3 |±1>=±|±1), (12) 

and the associated orthogonal projection operators 

P± = !(/ 2 ±<7 3 ), P±|±l) = |±l), P ± \ T 1} = 0. (13) 

Now we can rewrite as 

V(s) = E + ®P + U{s) + E_®P_U{s), (14) 

where 

U(s) = exp{i[sfiB A/2]a 2 } 

= cos(fj,B A/2)I 2 + issm(fiB A/2)(T 2 . (15) 



In (14) we have introduced <8> to separate explicitly the spatial and spin 



degrees of freedom. Adopting the standard matrix form for a 2 , we have from 



(15) the matrix representation 

U(s)=Af 



se 

-se 1 



(16) 



where 



e = tan(//B A/2) , N = 1/Vl + e 2 . (17) 



Here we have used cos(/ii?oA/2) > 0, which follows from the second of con- 
ditions (|8l). 



At this point we recognize (14) as the evolution operator for a QW on 
a line with unitary noise [HJ [9J [121 H3]- In that context, Hw is the 'walker' 
space and He the 'coin' space, E± governs steps by the walker to right or left 
on the line, and U(s) is the 'reshuffling matrix'. The walk is an essentially 
trivial one, with U = I 2 , contaminated by unitary noise that is characterized 
by the parameter e and the random variable s. 



6 



4 Analysis 

Depending on the direction of the magnetic field, the initial state vector 
10(0)) evolves to either 

|0(A)} = V(+l) |V(0)> or |V(A)} = V(-l) |V(0)) , (18) 

at time t = A, each possibility occurring with probability 1/2. But this 
simply means that the initial density operator p(0) = |-0(O)) (-0(0)1 = po, say, 
evolves into 

p(A) = Pl = § V(s)p V( S y . (19) 

s=±l 

More generally, at time t = NA, the density operator is 
p(iVA) = p N 

= ^n E V(s N )...V(s 2 )V(s 1 )p V(s i yV(s 2 y...V(s N )^ (20) 

Sl,S2, ... ,s N 

each s n being summed over the values ±1. 
Consider the case with 



|^(0)) = \k = 0) ® |r = +1) 



(21) 



so that 



Po = |0)(0|® 



1 









(22) 



using again the matrix representation as in (16). We postpone to the end of 
the paper a discussion of the difficulty of finding a positive energy state in 
this form for the system with Hamiltonian In this case, 



n+i)iv>(o)> 



NlE 



1 





e 






— e 



x |0) <g> 



1 





M(\l) 



1 




1) 



(23) 
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It follows that 



i/(+i) Po v(+i)t = W|i}<i|® 



-e 




1> < — 1 1 ® 



1 



-e 





(24) 



The expression for V(—l) po V(—ly is similar, with e replaced by — e through- 



out, and it then follows from (19) that in this case 



Pi=A/-(e) 2 |1><1|® 



1 






e 2 



(25) 



We see that pi, like po, is diagonal in the space of states spanned by all 
the vectors \k) <g> |r). It is easily shown by induction that the same is true of 
Pn for each non-negative integer N, and that in fact pn has the form 



N 
k=-N 



a Nk 

p Nk 



(26) 



where the prime indicates that the sum is over the N + 1 values — N , —N + 
2 , . . . ,N of k, and the constants c^at*:, are non-negative. 



It is also easily shown (T2J [13] from (20) that aNk, f^Nk are determined 
by the coupled recurrence relations 

a N+lk = N" 2 (aNk-! + e 2 f3 Nk -i) , k = -N + 1 ,-N + 3 , . . . ,N + 1 , 



(3 N+lk =U 2 {(3 N k + i + t 2 a N k + i), k = -N-l,-N + l,... ,N-1 



for N — 0,1,2 
with the initial conditions 

ot o = l, /3oo = 0, a fc = /3ofc = 0, fc 7^ . 



For example, (25) shows that 

a 11 =Af 2 , cn_i = 0, /3n = 0, /3 1 _ 1 =jV r 



•2 2 



(27) 
(28) 
(29) 
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Note that 



(+l\®{k\p N \k)®\ + l) = a Nk , (-l\®{k\p N \k)®\-l) = p Nk , (30) 

so that ctNk (resp. f3 N k) is the probability of finding the particle in the state 
\k) <8> | + 1) (resp. \k) <g) | — 1)) on measurement. Then 



P 



Nk 



Ct-Nk + @Nk 



(31) 



is the probability that the particle will be found in the state \k), with the 
value of t immaterial. Because the state \k) can be localized as closely as we 
like about x = kcA, we may say that P^ k is the probability of finding the 
particle 'at' that place, at time t = NA. 

When e = 0, the solution of (27) and (28) is onn = 1, with all other 
otNk and all P^k vanishing. The massless particle is free, and marches to the 
right at speed c, with probability Pnn = 1 of being at x = NcA = ct at time 
t = NA. If we had chosen r = —1 in the initial state, the particle would 
have marched to the left. In each case the motion is always purely ballistic; 
there is no transition to diffusional evolution at large times, because there 
is no noise contaminating the QW, which has the trivial reshuffling matrix 
U — li- If we were to choose an initial mixed state instead of (22), taking 
Po — |0)(0| ® ^2/2, the resulting probability distribution on the x-axis would 
have Pnn = Pn-n = 1/2, giving an extreme example of the familiar two- 
horned distributions associated with simple QWs on the line [T5] . 

When e = 1, p7fc and fl3ll) show that 



P 



N+lk 



Nk-l 



P 



Nk+1, 



(32) 



which is the defining rule for the Pascal's triangle of successive probability 
distributions centered on k = for a simple classical random walk (CRW), 
leading to 



Pnh 

for N = 0,1,2. 



2 N ^(N+k)/2 ) 



k = -N , —N + 2 . 



N 



(33) 



where C v = p\/q\(p — q)\ is the binomial coefficient. 



Then ([27]) and (J28J) give 

OlNk 



— r N ~ x k 

2 N ^(N+k-2)/2 ' ^ 



-iV + 2 , -iV + 4 , . . . iV . 



P. 



Nk 



iP> 



N-lk+1 



1 ynN-l 
2N Ly (N+k+2)/2 



k = -N ,-N + 2,...N -2 



= , P 



N N 



«0 



Poo = 0. 



(34) 
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Thus there is no ballistic regime in this case; the evolution is purely dif- 



fusional. As is well known [16], the distribution (33) is asymptotic to a 



continuous, normal one centered on k = as N — > oo, with density 

-k 2 /2N 

P(N, k) = , , -oo < k < oo . (35) 
V / 2A% 

Then we can say that after a time t = NA with N ^> 1, there is 95% 
probability that the particle can be found within two standard deviations of 
the origin, that is to say, with \x\ < 2a//VcA = 2c\/tA, in sharp contrast to 
the behavior of the free, massless particle, which would always be found at 
a distance ct from the origin. 

In terms of the physical variables of the relativistic particle, we have 
asymptotically as t = NA — > oo, the density 

e -x 2 /4Dt 

V(t,x)= . -oo<x<oo, D = c 2 A/2. (36) 



Note however that (|8j) requires e < 1, so that our approximate treatment of 
the relativistic system may break down when e = 1. 



For general values of e, the solution of (27) and (28) has been found in 
the form IT31 



a NN = Q(N,N), p NN = 0, 

and 

a Nk = Q(N } k) - Q(N -l,k+ 1)/(1 + e 2 ) , 

f3 Nk = e 2 Q(N — l,k + 1)/(1 + e 2 ) , (37) 
for k = —N , -N + 2,... ,N-2, where 

(N-\k\)/2 

Q(N,k) = — — CiN-Uy^^-iy. (38) 

^ ' s=0 



Then (31 ) gives 



Pnn — Q(N, N) , and P Nk = Q(N, k) - jj-^ Q(N - 1, k + 1) 

for jfe = -JV,-iV + 2,...,JV-2. (39) 
As N = 1 1 A increases for any given < e < 1 , the probability distribution 



(39) undergoes the transition from ballistic evolution to diffusive evolution 
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mentioned in the Introduction. This is seen in Fig. 1, which shows probability 
versus x/cA = k in the case e = 0.2, for JV = t/A = 20 , 50 , 100 , 200 . The 
plots reveal also that during the transition, the probability distribution may 
be considered to consist of two components: (1) the probability at the point 
x = ct corresponding to ballistic motion of the massless particle at speed 
c, marked by the circled point on the right in each plot, which decreases 
as time t = NA increases, and (2) a growing diffusional distribution on x- 
values closer to the origin, eventually swamping the ballistic component, and 
approaching a normal distribution centered on a mean positive displacement. 



This transition is reflected in the behavior of Pnn, the probability that 
the particle is at x = ct at time t = NA (the circled dot on the plots in Fig. 
1). This probability equals 1 in the free particle case e = 0. From (39) and 



(38) we have in general 



P 



NN 



,2\N 



[1/(1 + eX 3? 3? = V(l + O 



(40) 



As iV = t/A increases, there is a steady decline in the probability that the 
particle continues to travel at speed c to the right. This is clear in the 
successive plots of Fig. 1, and is shown explicitly in Fig. 2. 

It is also instructive to consider the moments of the probability distribu- 



tion (39), defined for each N as 



N 



S ( n )= ' P Nkk p , P = 0,l,2 



(41) 



-N 



The zeroth, first and second moments have been calculated exactly [13], as 



S 



(o) 

N 



s 



(1) 

N 



2e 2 



2\ N ' 



c(2) _ 1 



2Ne 2 - 1 + e 4 + 



'1 -e 



2\N+1 



(1 + e 2 ) N ~\ 
In the ballistic regime, iVe 2 1, and 
(1 _ e 2 )^+7(l + e 2 )^ 1 = 1 - e 4 - 2Ne 2 + 2iV 2 e 4 + 0([iVe 2 ] 3 ) 
so that 

S^ = ±(Ner + 0([Ner). 



(42) 



(43) 



(44) 
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The rate of growth of the second moment, quadratic in N, is characteristic 
of a QW [IS]. 

In the diffusive regime, Ne 2 3> 1, and 

= \ (Ne 2 ) - ^ + 0(l/[iV e 2 ]) , (45) 

showing the rate of growth is now linear in N, as typical of a CRW [16J. 
Note that however small is e 2 , eventually N becomes so large that the second 
regime is reached. The changeover occurs in the region where Ne 2 = 0(1). 
This is the behavior observed in numerical simulations of quantum walks with 
unitary noise [H], of which the present process may be considered a simple, 
special case, with the advantage that it is more amenable to mathematical 
analysis. 

Note also that the second moment can be interpreted as the expecta- 
tion value ((x/cA) 2 ). Fig. 3 shows the transition from quadratic to linear 
behavior of this quantity with increasing N = t/A in the case e = 0.2. 

The asymptotic behavior as iV — > oo of the moments of the distribution 



(39), defined as in (41), has also been calculated [13j, to give 



(2n) iV" (2n-l)! 
N ~ (n-l)!2"- 1 e 2 " ' 

(2n+1) N n (2n + l)! 1 - e 2 
Sn ~ "^!2^ W ' (46) 

for n = 0,1,2,.... From this it follows that as iV — > oo, the probability 
distribution given by P^k is asymptotic to a continuous, normal distribution 
with density 

e 2 (A ._ fco) 2 /2A r 

P(N,k) = . -oo<A;<oo, (47) 



where ko = (1 — e 2 )/2e 2 . This generalizes the result (35), which is recovered 
when e = 1. Fig. 4 shows the degree of agreement between the exact and 
asymptotic probability densities for e = 0.2 with N = t/A = 100 and N = 
300; in the latter case, the plots sit almost one on top of the other. Note 
that because the P^k are defined at A;- values two units apart, we have for 
each value of N, 



E 



\P Nk A k = 1 , A k = 2 , (48) 
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showing that it is the set of values k = —N , —N + 2 , . . . ,N that is 

to be considered a probability density for comparisons as in Fig. 4 with the 
continuous density P(N, k). 

For a general value of e, we now have that as t = NA — > oo, the proba- 
bility density on the x-axis is asymptotic to that for the normal distribution 

e -(x-x Q ) 2 /(4Dt) 

V(t, x) = , , — OO < X < oo , 

V4:irDt 

x = (l- e 2 )cA/2e 2 , D = c 2 A/2e 2 . (49) 
Here the displacement of the mean position (x) of the particle at large times, 



to the value Xq, is noteworthy. Using (17), we see that this value is given in 



terms of the physical variables defining the Hamiltonian as 

_ cos(/x£? A) 
X ° ~ 1 - cos(^oA) CA • (50) 

It can take any positive value, and is independent of the (very large value of 
the) time t. We can also say that as t — )■ oo, there is a 95% probability of 
finding the particle within two standard deviations of the mean; this is the 
probability to have \x — xq\ < 2c\/tA/ tan(fxB A/2). 



5 Discussion 

The transition from ballistic to diffusional behavior in QWs with unitary 
noise is surprising; in the context of a relativistic quantum system as we have 
discussed here, it is even more remarkable. Diffusion commonly arises as a 



non-relativistic, classical process. Note from (17) and (49) the dependence of 
the associated diffusion coefficient D on the speed of light as well as on the 
magnetic moment, the strength of the magnetic field, and the time interval 
between changes in field direction. 

Can a system of the type we have described be realized physically? There 
are several difficulties standing in the way. In the first place, it is now 
thought that all neutrinos have small rest masses, leaving us with no can- 
didate massless Dirac particles. Even if it should transpire that one of the 
neutrino types is after all massless, there is no evidence for neutrino mag- 
netic moments, although there has long been speculation that such might 
exist [HI HH HE 1201 EI] • 

Could the model apply to neutrinos, or indeed to the neutron, in situa- 
tions where the rest-mass contribution mc 2 (3 now properly appearing as an 
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extra term in the Hamiltonian (HI), can nevertheless be neglected relative to 
the magnetic interaction term, so that (Jll) still applies? Evidently, this would 
require 

\HfiB /2\ > mc 2 (51) 



in addition to the conditions (|8j). Inserting the observed values of \i and m 
for the neutron gives values of \Bq\ several orders of magnitude greater than 
any that have been achieved in the laboratory, although it is conceivable that 
such field strengths could occur in extreme cosmological situations. Note that 
according to (pi), extremely large \B \ values require extremely small field- 



oscillation times A if our approximate treatment is to be valid. In contrast 
to the situation for the free Dirac equation [7J, where the size of A can be 
adjusted at will to ensure accuracy of the approximations used there, in the 
present case it is determined once the external magnetic field is prescribed. 
For a neutrino (with anomalous magnetic moment) and very small rest 



mass, (51) is more easily satisfied. Although still unlikely in the laboratory, 
a process of the type we have described might apply in supernovas [22] or 
in the solar corona, for example, and affect significantly the fluxes of one 
or more neutrino types through the very strong stochastic magnetic fields 
occurring there [T71 [T8l [T9l [23] . 



A more subtle difficulty that has been raised following (22) concerns the 
form of positive-energy states of the relativistic particle. Just as for a free 
electron [21] , it is possible to construct positive-energy states of a free neu- 
trino that are arbitrarily highly localized about any given point [25J , and we 



could have used such a state in place of the \ip(0)) in (21 ), supposing that the 
particle is free at t — 0_, and that the interaction is switched on at t = + . 
However, such a state and its translates do not have compact support and 
are not mutually orthogonal, severely complicating the analysis of the QW. 
Such a more complicated analysis is not warranted, in our opinion; for even 
if the particle is in a positive-energy state at t = 0_, it will not be so at 
t = + when the field is switched on, because the Hamiltonian, and its pos- 
itive energy states, change form. Similarly, even if it can be arranged that 
the particle is in a positive-energy state at t = NA_, it will not be so at 
t = NA + if the field abruptly changes direction at that time, for the same 
reason. This difficulty, which is perhaps related to the Klein paradox [26], is 
not peculiar to the system we have discussed here, nor to our way of treating 
it. It seems clear that it must beset the analysis of a massless particle in any 
time-dependent, discontinuously changing external field. One could try to 
overcome the difficulty by projecting onto positive energy states immediately 
after each discontinuity occurs in the Hamiltonian, but it is not clear how 
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to do this in a simple way that preserves the unitarity of the resulting time 
evolution and hence the length of the state vector. 

We have not attempted to resolve this difficulty here, contenting ourselves 
with analyzing the model as described, in the belief that it makes an inter- 
esting addition to the collection of relativistic quantum systems that have 
been considered previously, and suggests that a quite new type of behavior 
can appear in such systems when subject to stochastically varying external 
fields. We hope that it provokes further study. 

It would be interesting, and more realistic, to consider the Hamiltonian 
(IT]) for the neutron, with rest-mass term added, and evaluate the evolu- 
tion numerically, without making any approximations like ([7]), to see if the 
stochastic nature of the interaction term continues to lead to diffusional be- 
havior at large times. However, the second difficulty mentioned above would 
still require resolution. 

The exact evolution in the massless case, associated with (|6]), could also be 
treated numerically to indicate any limitations of our approximate treatment. 

Finally, it should be mentioned that the present model can also be con- 
sidered in the context of quantum simulations of relativistic effects using 
the experimental apparatus of trapped ions. Indeed simulations of Dirac's 
equation and associated relativistic quantum effects for a single trapped ion 
have already been proposed [27] and experimentally realized [28], especially 
concerning the simulation of the Zitterbewegung phenomenon and the Klein 
paradox [291 US]- Simulations of a Dirac particle in a magnetic potential 
and its topological properties, using trapped ions, have also been proposed 
|31j . It is conceivable that the machinery of trapped ions might also allow 
a simulation of the Dirac QW driven by fluctuating magnetic fields that we 
have developed above, although it would be challenging for such simulations 
to overcome the physical obstacles indicated. 

Acknowledgment: We thank the School of Mathematics and Physics, Univer- 
sity of Queensland, for its hospitality during visits by D.E. (sabbatical) and 
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Figure 1: Clockwise from Top Left: Plots of probability against x/cA, for 
t = 20A,50A,100A,200A, with e = 0.2. The probability at x = ct is 
marked with a circled dot in each plot. Note the different scales on the axes 
in the plots. 
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Figure 2: Probability that the particle is at x = ct after time t, in the case 
e = 0.2. 
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Figure 3: Change from quadratic to linear growth of ((x/cA) 2 ) with increas- 
ing t = NA, in the case e = 0.2. The transition occurs in the neighborhood 
of N = 25, where Ne 2 = 1. 
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Figure 4: Comparison of the discrete probability density \PNk, k = 
—N , —N + 2 , . . . ,N, marked with dots, with the asymptotic normal form 
P(N, k), marked with circles, for e = 0.2 and iV = 100 , 300. Note the value 
at x — ct marked with an asterisk in the first plot. 
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